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Abstract 

This is intended to be a very basic introduction and short (and of course incomplete) 
overview of thermal field theory. In the first part, I introduce the thermal propagator 
at a very simple level and give the Feynman rules using the time-path contour method. 
In the second part, I give examples of these rules in scalar theory and discuss the 
origin of the thermal mass and other important effects as infrared divergences and 
phase transitions. In the third part, I outline the resummation program of Braaten 
and Pisarski. 
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1 Introduction 



The goal of Thermal Field Theory (TFT) is to describe a large ensemble of multi- 
interacting particles (including possible non-Abelian gauge interactions) in a thermo- 
dynamical environment [jl], |] . The main activity of physicists in that field deals with 
the description of the canonical or grand-canonical ensembles, that is of systems at 
fixed temperature and chemical potential 0. However, extensions to non-equilibrium 
systems also exist [4-7] (as well as a description in the micro-canonical ensemble [^), 
but I will not discuss them here. 

The TFT approach is different from older and more familiar theories such as the 
kinetic theory ^] or many-body theory |11| in the sense that the TFT is endowed 



with the following advantages: i) it uses a path integral approach, ii) it can treat 
non-Abelian gauge interactions as QCD and iii) it is Lorentz-covariant. 

Although the TFT is now almost 40 years old, its usefulness has been particularly 
acknowledged in the middle of the 70s, for the study of phase transitions in quantum 



field theories. Lately, most of the developments of the TFT [|I2], |T3| have been moti- 
vated by the study of QCD at finite temperature, as the creation of the quark-gluon 
plasma in the laboratory becomes feasible. 

Within the TFT, two classes of formalisms can be distinguished: one class is based 



upon a complex-time contour 0, ^ |T^, |T^, |T6] and the other on C* algebra . I will 
only describe the former. 

The main idea of the TFT is to use the path-integral approach of the usual vac- 
uum field theory and describe the temperature that appears in the Boltzmann factor 
exp(— /5iJ) using the time-evolution operator with complex time variables t = —i(3 = 
—i/T (everywhere conventional units h = c = k = 1 are used). 



Historically |T8[, Matsubara |T9| was the first to build a TFT by incorporating a 
purely imaginary time variable into the evolution operator. His name is associated 
with the discrete energy frequencies, the so-called Matsubara frequencies onto which 
one has to sum over in the Imaginary-Time Formalism (ITF). Then came some im- 
portant contributions from Schwinger |T^, Mills and Keldysh |^ who developed a 



formalism based upon the choice of a contour in the complex plane. This is called the 
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Real-Time Formalism (RTF). The latest developments include a functional formula- 
tion of the theory [|l^. Independently, Umezawa and coworkers ^ have taken 
a different approach, based on C* algebra, called Thermo-Field Dynamics (TFD), 
which turns out to give essentially the same results. 

The main applications of the TFT in high energy physics can be sorted in three 
classes: 

i) Cosmology: Obviously, if one is interested in the studies of hot plasmas, the 
early universe is a very good example. At any time before the recombination, the 
expansion length of the universe is much smaller than the mean free path between the 
particles, so that one can really speak of a perfectly thermalized plasma (at least under 



certain conditions [|22[). This is also clearly demonstrated by COBE's results on the 
cosmic microwave background radiation, which shows a purely uniform black-body 
spectrum, up to fluctuations 6T/T ~ 10^^ p3[. The temperature of the primordial 



plasma ranges from the Planck mass down to a few eV p4| . 

In fact, the first and major success of the TFT concerns the symmetry restoration 
in spontaneously broken gauge theories at finite temperature [25-27]. This is also at 
the basis of the inflationary scenario [^]. This topic has recently been revived in the 
case of the electroweak theory, where the order of the phase transition is crucial for 



baryogenesis [29 



There exist other types of study in the early universe, as for instance the calculation 
of reaction rates taking place in the hot plasma in order to determine the abundances 



of some species |^0[ . One should also mention an excellent work by Rebhan on thermal 



corrections to quantum gravity |31 



ii) Astrophysics: Cores of neutron stars, supernovae, red giants and white dwarfs 
are composed of extremely dense plasmas (p = 10^-lO^^g/cm^). The temperature 
can also be very high during the collapse of a supernova (a few tens of MeV). There 
does not exist many applications of the TFT in such systems, but recent studies 
on neutrinos and axions emission rates in these stars have proved to be quite useful 
[32-38]. 

iii) Heavy-Ion collisions: This concerns the formation of the quark-gluon plasma 
(QGP) in the laboratory [|39|. Lattice calculations clearly show that QCD has a phase 



transition at Tc = 100-200 MeV, above which hadronic matter is deconfined [^]. The 
TFT is an obvious and unavoidable tool in the study of the QGP and a lot of work 
has already been done in the calculation of signals from such matter [41-43]. 

This is almost all concerning particle physics. Of course, there are other fields of 
application as for instance condensed matter, for which TFD was in fact developed. 

To be honest, the plasmon decay into vv pairs [^, which is the dominant cooling 
mechanism of red giants and white dwarfs, is so far the only experimental verification 
(although indirect) of the TFT in high-energy physics. Therefore, the support of such 
theories proceeds more by analogy with what is known from the classical case. It is 
hoped that each field cited above will bring, in the near future, observable verifications 
to these theories. 



2 Thermal propagator 

In this section, I will derive the thermal propagator in the simplest possible model, that 
is a free scalar field theory. All the relevant features of the TFT show up already at this 
simple level, so that extensions to more sophisticated theories will be straightforward. 

The starting point of the TFT is the new definition of an observable A for a system 
in contact with a heat bath at fixed temperature T (in the canonical ensemble): 

(2.1) 

where Z is the partition function. The Boltzmann factor weighs the occupation num- 
ber of the states that are accessible to the system (here the trace is performed onto 
discrete states, but it can of course be continuous). The next step is to realize that 
the Boltzmann factor can be described by a rotation of the operator in the Heisenberg 
representation 

e'^^A(t)e-''^ = v4(t-i/5), (2.2) 

where t is now a complex time variable. For analytic reasons, only some complex time 
values are allowed. For instance, by evaluating the double product 

Tr e-'^^A(t)B(t') = ^e-^^'"(m|A(t)5(t')|m) 

= ^e-'^"(*-*')e'^'"(*-*'+'^)(m|A(0)|n)(n|S(0)|m), (2.3) 
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it appears that the trace operation is meaninful only for time arguments in the range 

-(3< Im(t - t') < 0. (2.4) 

Consider now the 2-point Green function in a scalar theory. By definition 

G{X-Y) = (T^XmY)))^ 

= 9,{xo-yo)G+{X-Y) + e,{yo-xo)G_{X-Y), (2.5) 

where a complex time ordering has been defined {6c{xo — Vo) = 1 if l/o precedes xq 
along the contour) . The domain of analyticity of the 2-point Green function is given 

by 

- P < lm{xo -yo) <0 when 9c{xo - yo) = 1 

< Im(xo -yo) <P when Odyo - Xq) = 1. (2.6) 

This means that the contour in the complex plane must not go up but can only go 
down. 

An important property, the Kubo-Martin-Schwinger (KMS) relation, can also be 
derived: 

G+{xo - (yo + ^/3), ^ - = G^X - Y), (2.7) 
which follows from the trace invariance by circular permutation. 

We are now in a position to write down the 2-point Green function for a free 
scalar field in terms of operator kernels. Consider a superposition of creation and 
annihilation operators in the plane-wave approximation 

= / TV^^^^'^^y''''' + «(^)^"^^-'']' (2-8) 
J (27r)2 ^/2uJ 



where uj — ko — ^/\^\^~-\^m? . One finds for the thermal propagator 
G(X -Y) ^ f / (2 9) 



X 



I e,{xo - yo) [e-^-^+^^-^(a(X)at(P))^ + e^^-^-^-^(at(P)a(X)) 
+ 0c{yo - xo) [e^^-^-^^-^(a(/r)at(P))^ + e-^^-^+^^-^(at(i^)a(P));j] } 



There remains to calculate the thermal kernels. Taking for Hamiltonian 

: a{K)a\K) + a\K)a{K) :, (2.10) 

one can start from the thermal average of the commutation relation 

{a{P)a\K))p - {a\K)a{P))p = 5'ip-k). (2.11) 

And using the cyclicity of the trace and the commutation relation [H, a] = —ua, one 
arrives at 

{a\K)a{P))p = nB{uj)6\p-k) 

{a{P)a\K))p = {l+nB{u))5\p-k), (2.12) 
where I have defined the Bose-Einstein statistical weight 

^b{^) = ^p^Y (2-13) 
The propagator can then be rewritten as 



G{X -Y)= I e^{x, - yo)e-^^-(^-^) + e^{y, - x,)e^^<^-''^ 



(2.14) 



At this level, one has to specify the contour. As the contour must start from some 
initial time and go down to to ~ T'I^i the simplest choice is just a straight line along 
the imaginary-time axis. This choice of contour leads to the ITF and is the oldest 
and most widely used formalism in the TFT. Parametrizing xq — yo = one has 

G(r, x-y)= / ..^ e-^'-("-^^) (1 + nB(cu))e-"l"l + nB(cu)e"l"l . (2.15) 



In momentum space, the theory reduces to a 3-dimensional euclidean theory with 
an infinite summation over the Matsubara frequencies ujn = i27mT,n g] — oo,+oo[. 
Explicit calculations show however that is more convenient to stay in r-space . 



The obvious disadvantage of this contour is to lose completely the real-time argu- 
ment. In principle, with this formalism, one is restricted to calculate static thermody- 
namical quantities as the free energy. This has motivated several people to consider 
a different contour, which would keep a real-time argument. 
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The contour of the RTF is shown in fig. 1 and is composed of four different pieces. 
This is the simplest choice of contour if one wants to have real-time arguments. For 
analytic reasons, it can be shown that the contributions from C3 and C4 can be 
neglected One is therefore left with two possibilities for both xq and yo to lie 

either on Ci or C2. This gives four different propagators, which are usually written 
in a matrix. They are labelled as "11" propagator, corresponding to Xq € Ci and 
Ho G Ci, and so on. An appropriate choice of the parameter cr is the intermediate 
value a = (3/2. In this case, one finds the following form of the complete propagator 
in momentum space 



/ 



G^\K) G^\K) 



UiP,K) 



^ A{K) 



\ 







G{K) 
where 

/ 

sinh 9k cosh 9 k 




UiP,K), (2.16) 



^ cosh 9k sinh 9k 



u{p,K)= ;'77 ^"77 I, (2.17) 

and 

cosh^;^ = (1 - e-^")~'^' , sinh^^ = e"^"/' (l - e"'^")"'^' . (2.18) 
In (2.16), A is the Feynman propagator at zero temperature 

= r.2 \ ^ • • (2-19) 

j^z _ _|_ 

The welcome surprise is that these are just exactly the propagators in TFD. This is 
a very strong support of the theory in the sense that TFD and the RTF do not use 
at all the same mathematical framework. 

In TFD language, the 2-type propagators: the "12", "21' and "22", are associated 
to "ghost" fields (not to be mistaken for Fadeev-Popov ghosts!). They are unphysical 
since one of the time arguments has an imaginary component. The only physical 
propagator is the "11" component 

G^\K) = A{K) + 2TTnB{iu)6{K^ - n?). (2.20) 

This propagator was derived in the early attemps at real-time techniques It 
clearly shows the finite-temperature contribution in the (5-piece. But, in order to have 
a consistent theory, the other propagators must also be taken into account. 



6 



The Feynman rules in the RTF are the following: diagrams have the same topology 
as in the vacuum theory, and the same symmetry factors. Furthermore, there are two 
types of vertices, 1 and 2. Physical legs must always be attached to type-1 vertices, 
never to type-2 ones. One must sum over all the configurations of type-1 and type- 
2 vertices according to the above rule. Then, one has to use a G'^^{K) propagator 
between a vertex of type a and a vertex of type h when K flows from a to h. All the 
rest remains the same as in the vacuum theory. 

The Feynman rules shown above are those of a free-field theory. Implementing 
interactions is rather straightforward and is done in the same way as in vacuum, as all 
the necessary formalism is already set up. For cubic, ((?/3!)0^, or quadratic, ((7/4!)0'*, 
interactions, the Feynman rule is to attach —ig at a type-1 vertex and +ig at a type-2 
one. In fact, in order to construct these rules one has to make the assumption that the 
interaction is adiabatically swichted on and off around t = 0, which is just the usual 
assumption in vacuum theory [^. In a heat bath, this is certainly questionable, as 



one might think of any particle in constant interaction with its thermal surrounding. 
One can phrase the question differently by asking: are there asymptotic states in 
a heat bath? Indeed, practical calculations show some difficulties and formal 
considerations seem to give a negative answer Still, perturbative calculations 



can make sense, as I will try to show in the next sections. 

The field theory at finite temperature is renormalizable, provided the vacuum 
theory is so. This is intuitively obvious as the thermal corrections all come with a 
Boltzmann factor that cuts off any ultraviolet divergence. At the perturbative level, 
the problem appears more subtle, but is clearly not an issue. 



The path-integral formulation of the same theory is a straightforward exercise [IB 
The generating functional is given by 

Z\j] = Z[0](T,exp (t j^d^xj{x)(t>{x)^)p. (2.21) 

and the path is such that x{t — i(3) = x{t). By differentiation with respect to the 
sources, one obtains the n-point Green functions: 

1 6^Z[j] 



G (Xi , . . . , Xn) 



Z[0] i6j{xi) . ..i6j{xn 
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(2.22) 

j=0 



With the contour C shown in fig. 1, the same matrix propagator as in eq. (2.16) is 
easily recovered. For completeness, I also give the gauge boson "11" propagator in 
Feynman gauge 

Gf,{K) = -g,, (AiK) + 27rnBico)5iK')) . (2.23) 

Similar relations hold for the other components (one just has to multiply the thermal 
matrix (2.16) by the metric). For Dirac fermions, one has 

S'p\K) = i^^K^ + m) [A{K) - 2nNF{ko)S{K'' - m^)) , (2.24) 

where the Fermi-Dirac statistical weight has been introduced (in presence of a chemical 
potential fi) 

^P^''^> - e/3(fco-M) + 1 + e/3(-fco+M) + i- ^^■^^> 
Again, one has to multiply the Dirac structure to the fermionic scalar matrix in order 
to have the other components. One should also notice that the ghost gauge fields, 
although anticommutant, obey the Bose- Einstein statistics. 

After some initial confusion, there has been a considerable amount of work in the 
past few years in comparing the different formalisms. The first remark is that, by using 
the RTF with a free parameter a, one gets the TFD for a = (3/2, and the Keldysh, or 
also called the "cut" propagators, for a = or cr = /5. The comparisons between RTF 



and ITF are much more complicated, since the Feynman rules are rather different 
However, there is now a general agreement that they give the same physical answers 
3^]. In fact, it is also possible to calculate dynamical quantities in the ITF, by analytic 



continuation to real-time values uj ^ ko + ie. Note that the equivalence between the 
two formalisms is obtained by taking the analytic continuation uj —>■ kQ + iek^. The 
question of which formalism is the easiest to use is a difficult one. Each has its 
advantages and its weaknesses, and using one or the other is probably more a matter 
of taste than of convenience. However, one should note that when dealing with non- 
equilibrium systems, a real-time approach is mandatory. For systems at equilibrium, 
the best solution is certainly to use both formalisms, as this provides a useful check 
on the calculations. 
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3 Some examples in scalar theory 



3.1 Tadpole diagram and comparison between ITF and RTF 

The simplest thing to calculate in the TFT is the tadpole diagram in a scalar the- 
ory with quadratic interactions as gcj)^ theory. Is may also be the most instructive 
calculation. The Feyman rules in the ITF give 

g ^ [ d?k 1 

where uJn = 2imTT. In order to perform the summation over n, the standard trick is 
to use a contour integral 

T 2 /(fco = ^u;^) = ^ / dA;o/(fco)^ coth (3.2) 

n=-oo II 

provided /(A^o) has no poles on the imaginary axis. The contour can be deformed in 
such a way that (see fig. 2) 

+2P 1 r r+ioo r — ioo 1 

T f{h = lUn) = — / dkof{ko)+ dko{f{ko) + f{-ko))^^—- . 

(3.3) 

This allows us to separate the vacuum part (first term) from the finite-temperature 
contribution (second term). Using the explicit value of /(^o) = ~ ^o) with 

uj = vP + m^, one gets for the thermal contribution to the tadpole diagram 

g f d^k 1 



This result can be derived in a much simpler way by using the "mixed" propagator 

m 

„ +00 j3l „/3 

^ ^ "^'^ ^^dr5(r)A(r,u;) 

l + 2nB{uj)). (3.5) 



2 J (27r) 
g f d^k 1 



2 J (27r)3 2cj 

The "1" piece is again the vacuum contribution, which is quadratically ultraviolet- 
divergent while the "272^" piece is the finite-temperature contribution. 
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With the RTF, the separation from T = to T 7^ is automatic from the begin- 
ning, so that one can write 

= f / ^^^'^^^^ " m')ns{u;). (3.6) 

The integration on ko can be trivially performed with the help of the 5-function and 
the same result is obtained as with the ITF. 

This one-loop example does not illustrate well the complexity and the relative 
advantages of each technique, which are observed at higher orders. One can just see 
that the methods are quite different. 



3.2 High-temperature expansion and thermal mass 

With this one result in hand, let me continue the discussion. We have 

The remaining integral cannot be computed analytically, except in some limiting cases. 
Let me discuss the high-temperature limit, which is of greatest interest. One has 



The tadpole diagram is the first-order correction to the free propagator and there- 
fore the above result is nothing but a mass term (because it is independent of the 
momentum). In the limit of zero bare mass, one has the contribution 

m} = (3.9) 

This is an effective mass and corresponds to the fact that in the heat bath, the 
propagation of particles is altered by their continuous interactions with the medium. 
This is certainly one of the most important results of the TFT. 

The fact that even massless particles acquire a mass does not affect the chiral 
symmetry. The thermal mass is generated radiatively and comes from a term that is 
chiral-symmetric if the initial Lagrangian is. 
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3.3 Higher orders and infrared problems 

In order to see if the perturbation series is well behaved, it is necessary to look at 
higher orders. At two- loop order, two different topologies contribute (see fig. 3). 
According to the RTF rules, one has for the first topology 

^^a\K) = / ^,Dl,iP) J ^ADUK) - (3.10) 

where I have used the fact that the type-2 tadpole is identical to the type-1 tadpole. 
It is interesting to note that 

D^K) - DUK) = A\K) - 2PP / ^ 2n6{K' - m')nB{u;). (3.11) 

Thus, even if there are multiple products of statistical weights at intermediate stages 
of the calculation, they disappear when combining the different type-1 and 2 vertices. 
Also disappear the ill-defined 5-products. This is a generic feature of the RTF. In fact, 
a general theorem can be stated: there appears at most only one statistical factor (or 
a sum of) in a loop integral. This is obvious in the ITF, because it is the summation 
over the Matsubara frequencies that gives the statistical weight. 

As usual, the ultra-violet singularities in eq. (3.10) can be eliminated by choosing 
some suitable counter-terms (in fact the same counter-terms as at T = 0). Notice, 
however, that all topologies must be added in order to cancel the ultra-violet diver- 



gences, with a non-trivial mixing between T-dependent and T-independent terms 
In the high-temperature limit, the contribution to the thermal mass becomes 



6ml'^'^ = -^T' (- + ...), (3.12) 



and is infrared-singular when m 0. This power-like singularity originates solely 
from the first topology in fig. 3. Although infrared singularities are inherent in al- 
most all perturbation theories, whether at zero or finite temperature, there are good 
reasons to worry about this one: i) the above contribution enters as a correction 
to the thermal mass, and therefore participates into the solution of the pole of the 
propagator, which is supposed to be a physical quantity [Q; ii) at zero temperature, 
famous theorems as the Kinoshita-Lee-Nauenberg theorem Q demonstrate that 
singularities appearing at intermediate stages of the calculation cancel in the final 
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physical result. However, a crucial argument of the Lee-Nauenberg version of the 
theorem is the finiteness of the mass shift at all orders ||5^, which is just the quantity 



considered here; iii) indeed, this mass shift will enter into physical quantities as it acts 
as a kinematical cut-off and will appear, say, in the production rate of light weakly 
coupled particles from the heat bath. 

So, what can be done with this infrared singularity? If it cannot be cancelled, then 
one solution is to use an infrared regulator. In fact, an obvious infrared regulator is 
just the thermal mass obtained at first order, m^. In that case, the resulting behaviour 
is 

5ml^'^ = 0{g^T^). (3.13) 

That is, the result is still perturbatively acceptable as it is smaller than the previous 
correction. Unfortunately, this hope does not survive higher-order corrections. Indeed, 
by attaching an arbitrary number of tadpoles to the first line (see fig. 4), one gets 
always the same correction 

Sml^""^ = 0{g^gT'). (3.14) 

An infinite number of diagrams contribute to the same order in the coupling constant. 
This clearly shows the breakdown of perturbation series at finite temperature. 

On the other hand, this infinite subset of the most infrared diagrams can be 
resummed. It is nothing but a Taylor expansion 



mf ( d kdu 



oo 

^0 N\ \drr?j Jm e^^ 



°° k du 

m=0 Jmp — . 



(3.15) 



The right-hand side of this equation is perfectly defined for any value of the coupling 
constant g. Seen from this perspective, the infrared problem originates from the bad 
expansion of y/g around zero. 

The resummed result is not very different from the first-order correction, m^, 
which, in the light of the above equation, is not surprising. This seems to make 
perturbation theory at finite temperature reliable. I will come back to this important 
point in the last section. 
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3.4 Phase transitions 



The appearance of a thermal mass has another important consequence. Consider a 
bare Lagrangian given by 



m = \{d,<Pf + \li'<P'-^<p\ (3.16) 



This Lagrangian has a negative square-mass term. The minimum of the effective po- 
tential is degenerate and occurs for non-zero values of the field. Now, the temperature 
introduces a positive square-mass term. The effective potential at high temperature 
then looks very different 

as it has only one minimum at 0c = 0. Hence, there must be a phase transition 
between the high- and the low-temperature regime. Whether this phase transition 
is of first or of second order is another question, which requires to take into account 
some sub leading terms Still, this is an important result, which has truly modified 



our way of thinking about the early universe. Any symmetry that is spontaneously 
broken in our present world, which is a low-temperature system (T = 3 K), will be 
eventually restored if going to sufficiently early times. This is true for instance for the 
Higgs mechanism, which is responsible for all the particle masses. Above T > 1 TeV, 
all particles become massless. 

This feature is not specific to gip'^ theory. With gcjr' interactions in a 6-dimensional 
space-time, the theory becomes completely unstable above a certain critical temper- 
ature, which can be computed exactly Phase transitions do also exist in gauge 



theories 27, 29 . 



4 The resummation of Braaten and Pisarski 
4.1 The idea 

We have seen in the previous section that in order to get rid of infrared problems, one 
has to resum a certain class of diagrams. However, as shown above in the case of the 
thermal mass, it is not always necessary to perform this resummation if one is not 
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interested in higher-order corrections. In the following, I will present a set of simple 
rules that determine in which case one has to resum infrared diagrams. The method 



is due to Braaten and Pisarski |5^, |53] and works also for gauge theories. 

Consider g'^cj)^ theory at very high temperature, so high that all other scales as for 
instance the bare mass, can be neglected. Still, possible other scales can be generated 
radiatively as gT, or "non-perturbatively" as {\n.g)T. The starting point of the BPR 
is to make a distinction between hard momenta of order T and soft momenta of order 
gT. The second step is to realize that the thermal mass is generated by a loop integral, 
where the momentum running inside the loop is hard. This is clearly what gives the 
(7^T^ contribution to m^. Soft momenta give higher-order contributions as g^T"^ terms. 
The hard momentum contribution to the loop integral is called a hard thermal loop 
(HTL). The third and final step consists in resumming only soft lines. Indeed, the 
corrections to the bare propagator being of order (7^T^, they start to be relevant 
only when K ~ gT, which is the soft scale. Hard lines do not need to be resummed 
and one can still use the bare perturbation series. 

With these simple rules, one has an improved perturbation theory. For an n-point 
Green function, the dominant contribution always comes from the HTL, which is of 
order g^T"^. If there is no HTL inside the loop, the integral must be cut into two pieces, 
a first piece where the momentum running inside the loop is hard, for which one uses 
bare propagators, and a second piece where one has to use resummed propagators 
as the loop-momentum becomes soft. The final result should be independent of the 
arbitrary intermediate cut-off of order y^T, which is put by hand to separate the two 
pieces. 

In g'^cj)'^ theory, there is only one HTL, which is the tadpole diagram. Therefore 
only the 2-point Green function needs to be resummed. The effective Lagrangian that 



leads to the improved perturbation series can be written as |^ 



i^eff = \{d,ct>f-\nil4>\ (4.1) 

with given by eq. (3.9). 

The very nice surprise of the BPR is that these features survive when considering 
gauge field theories. In QED, the 2-, 3- and 4-point Green functions have an HTL. 
The same is true in QCD, where the n-point Green function with only external gluon 
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lines, or with n — 2 gluon lines and 2 quark lines, must also be taken into account. In 
these cases though, HTLs are not just scalars but complicated analytic functions of 
the external momenta, g"'T'^f{uJi, ki). Another surprise is also that these functions are 
the same for QED and QCD (at least for the same number of bosonic and fermionic 
external lines). Finally, there are two extra bonuses: i) the HTLs are gauge- invariant 
and ii) they obey simple Ward identities. 

The HTLs can be derived from an effective Lagrangian given by |^ 



where () refers to an angular average; K is the hard loop momentum; the two param- 
eters uoq and are the plasmon frequency and the thermal fermion mass respectively 
(of order gT). With this Lagrangian, one has an improved perturbation series defined 
in terms of effective propagators and vertices (they are pictured in fig. 5 in the case 
of QED). 

In the next section, I shall present and discuss the HTLs for the bosonic and 
fermionic self-energies. 

4.2 The bosonic hard thermal self-energy 

An important consequence of the fact that the plasma constitutes a privileged rest 
frame is that the polarization operator can be decomposed into two different prop- 
agating modes, transverse and longitudinal (contrary to the vacuum case, both are 
physical): 

U^,{K) = P^MK) + Q,,IIl{K). (4.3) 
Expressions for the projectors can be found in The HTL is obtained by taking the 



high-temperature limit and by considering a soft external momentum, u,k <^T. After 
renormalization of the ultra-violet zero temperature piece, the thermal corrections for 
each mode are 

na.-)^3.2(i-|!)(:-^.^], (4.4, 
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where ujq is given by 

for an SU (N) gauge theory with Nj fermions. For QED, one has to substitute {Nf/2+ 
N)g'^ — > e^. Identical expressions are obtained in the uhra-degenerate hmit, fi ^ T, 
except for ujq, which scales with /i instead of T. When the fermion mass is not 
negligible, these expressions differ only slightly 

Of course, the above results are nothing new for QED, where they were derived 
already a long time ago using the kinetic theory . In that case, ujq is the plasmon 
frequency^ at which the electric charges oscillate in the plasma around their average 
positions. Because of this analogy, it is possible to call uoq the plasmon frequency for 
the general SU{N) gauge theory case. In fact, it is quite remarkable that expressions 
(4.4) are so general. 

In arbitrary covariant gauge with parameter ^, the resummed gauge propagator is 
(in the ITF) 

The two poles of this resummed propagator define two propagation modes (plasmons), 
which, because of the complicated structure of eqs. (4.4), are described by non-trivial 
dispersion relations. They are plotted in fig. 6. Notice that when k ^ 0, they both 
tend to the same value, u = ujq, which can then be regarded as a plasmon mass. 

There are also contributions below the light-cone, as the HTLs develop an imagi- 
nary part there. In the static limit, one has 

ReUTiuj^O,k) = kl — , 

ReULiuj ^0,k) = kl, (4.7) 

where k^y is the Debye mass corresponding to the screening of static electric fields 
(again following the same analogy with QED). Indeed, one finds that in non-degenerate 
plasmas kjj = e'^N^/T and in degenerate plasmas kjj = {e^ /■n'^)npF- For ultra- 
relativistic plasmas, = ScJq. Notice that fc/j also appears in the expression for 
the transverse propagator. However, it vanishes when — ^ 0. Hence, there is no 
screening for space-like transverse gauge bosons. This has a very important conse- 
quence as the resummation of Braaten and Pisarski is essentially based upon the 
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presence of a mass term of order g'^T'^ in the 2-point HTL. The transverse (or better 
called magnetic) mass can only be of higher order, that is m^ag = 0{g^T). In fact, it 
is exactly zero at all orders for QED |^. This is at the origin of the infrared problem 
in hot gauge theories l6l|, For an SU{N) gauge theory (A^ > 1) and beyond some 
order in the perturbative coupling constant, an infinite number of graphs contribute 
to the same order, due to diagrams with multiple interacting space-like transverse 
gauge bosons. In QED, or for a gauge theory at finite density, there is always enough 
screening coming from the fermionic loops so the problem does not seem to be so 
worrying. However, we shall see in the last section that there might be an alternative 
to this problem. 

Quite unfortunately, the literature contains may examples of confusion between 
the Debye and the plasmon masses. They are totally unrelated. The former acts as a 
screening mass for space-like longitudinal photons (for instance in a fermion scattering 
on a target). The latter damps the boson density in the plasma by a factor e~^°^'^. 



4.3 The fermionic hard thermal self-energy 

Following the same procedure as in the last section, the fermionic HTL is found to be 

^ f 1 , uj + k ( uj , uj + k\ k.7\ 

where mp is another thermal mass given by 

m'p = cp^, (4.9) 

and cp = {N"^ — 1)/{2N) is a color factor. In QED, one has to substitute cpg'^ e^. 
By looking at the pole of the resummed propagator, 

one finds, unexpectedly, two solutions. They are plotted in fig. 7. The first solution 
(+) corresponds obviously to the modified fermion propagation in the plasma, but 
there does not seem to be any intuitive physical interpretation for the second solu- 
tion (— ), apart from the fact that it arrives as a pure collective effect |^3|]. This is 
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a completely new feature, which was not seen before the TFT was used [^. Unfor- 
tunately, there seem to be no important phenomenological consequences of this new 
mode, sometimes called "plasmino" [ |65|1 . 



4.4 Damping rates 

The resummation program of Braaten and Pisarski was in fact developed to solve the 
long-standing problem of getting a gauge-independent answer for the gluon damping 
rate [|66l . 

The damping rates are computed through the discontinuity of self-energies at 
the pole of the propagator and have the following physical role: consider a particle 
distribution that is slightly out of equilibrium. One has 

dn{uj, t) 



dt 



-n{uj, t)Ta{uj) + (1 + an{uo, t))T,{uj), (4.11) 



where Ta and Fe are the absorption and creation rates of the given particle, respec- 
tively. Notice the different statistical factors for incoming and outgoing particles, with 
the parameter a to distinguish bosons (a = 1) from fermions (cr = —1). 
Equation (4.11) has for general solution 

n{uj, t) = + C(^)e-(^"-'^r^)*, (4.12) 

where C{uj) is an arbitrary function, which does not depend on time. Creation and 
absorption rates are related by the relation 

F,(a;) =e'^'^Fe(a;), (4.13) 

which is nothing but the KMS relation, written in a different form. Using (4.13) we 
have 

n{uj,t) = ^^ + C{u)e-^^\ (4.14) 

where 7 = (F^ — crFc)/2 is defined as the damping rate. Therefore, its physical 
interpretation is rather clear. It represents the inverse time scale it takes for a thermal 
distribution to reach equilibrium. In particular, the sign of the damping rate is of 
crucial importance as, if negative, the system is obviously unstable. 
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One can therefore imagine the controversy, which has lasted for many years, as 
the sign of the gluon damping rate was found to be dependent on the gauge choice! 
Braaten and Pisarski have shown that in order to get a consistent result, one has to 
use their effective perturbation theory . 



How to compute damping rates? The creation and absorption rates are related 
to the imaginary part of the self-energy [B^. In the RTF, Kobes and Semenoff have 



extended the cutting rules (or Cutkosky rules) to the finite temperature case ^8 
Again, results are similar in the ITF and in the RTF. 

From the results shown above, eqs. (4.4) and (4.8), it is easy to see that there is 
no HTL for the damping rates. Indeed, the HTLs have an imaginary part only below 
the light-cone, whereas the pole of the propagator is clearly above it. Therefore, the 
damping rates must be of higher order, g'^T, which makes them difficult to compute. 
Indeed, in order to go beyond the leading order in the BPR, one has to subdivide 
diagrams into different pieces, according to whether the internal lines carry hard or soft 
momenta. For soft external gluons, one has to use effective vertices and propagators 
(see fig. 8). For transverse and longitudinal gluons at rest (!) the result is 

7^ = 7^ - 0.264013/T. (4.15) 

So the sign is positive and independent of the gauge, which implies that the quark- 
gluon plasma is perturbatively stable. 

The case of a hard external gluon turns out to be very interesting, as there appears 
a singular behaviour that can only be cut by introducing a higher order cut-off of the 
type rrimag or 77- itself In this case, effective vertices are not used and only one 
soft internal line is resummed. Only the transverse gluon damping rate diverges and 
is found to be 

7T = ^^^ In^ -J + 1.09681... . 4.16 

Hence, jt appears both on the left-hand side and on the right-hand side of the equa- 
tion. In order to get this result, one has to use a resummed propagator also for the 
one hard internal line, which has nothing to do with the resummed propagator of 
Braaten and Pisarski. In fact, this propagator does not resum the HTL, but some 
other contribution. The result is also sensitive to the magnetic mass. This is bad 
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in the sense that nobody knows how to calculate this parameter consistently, except 
by going on the lattice. However, this uncertainty shows up only at the logarithmic 
level and has therefore very little phenomenological implications. The coefficient in 
front of the logarithm is on the other hand perfectly well known, and is consistently 
obtained by using the BPR. 

The damping rates for fermions have similar features. What of hot QED, where 
the magnetic mass vanishes? There, it seems that there is no consistent way to get 
the damping rate |6^, |70| and the problem is still open. Among the possible different 
solutions, the simplest could be that there is no pole at all in the fermion propagator 
[n]] ! It is also legitimate to ask the physical relevance of a quantity that is defined at 



a complex energy! On the other hand, one may note that the creation and absorption 
rates of fermions travelling through a heat bath are perfectly defined and are of order 
g'^T\ia{l/g), as expected W2 . 



4.5 Final remarks 

When computing at finite temperature, there appear singularities as in any other 
quantum field theory. Ultra-violet singularities are the same as in vacuum and are 
disposed of by renormalization. Infrared singularities occur when a momentum is 
vanishing and collinear divergences when a line goes on shell. According to the pre- 
vious discussion, the singularity can be screened by using a resummed propagator or 
vertex. However, terms that generate infrared singularities are not the same as those 
that generate the HTLs. Indeed, a typical on-shell expansion of the self-energy is (at 
first order in g) 

E(A-) = a/r= / + bi-K' I ^|-^, (4.17) 

where a and b are numerical factors. The ffist term is proportional to and rep- 
resents the HTL. The second term is infrared-singular and a priori infinite if one 
does not regularize in some way. However, by resummation, we know: i) that there 
are thermal cut-offs of order gT and ii) that the particle on-shellness is restricted to 
~ g^T"^ whether it is soft or hard. Therefore, the second contribution is of order 
(jf^ X g^T'^[l/g) ~ g^T"^ and is smaller than the HTL. It is clear then, for the BPR to 
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work, that the term with the screened singularity must contribute less than the term 
of the HTL. This happens to be true indeed, except of course in the case where there 
is no screening, i.e., when the soft line is a magnetic boson. On the other hand, it has 
also been found that those kinds of infrared singularities can cancel out in the physical 
answer. Indeed, although a complete general proof is still lacking, there is by now 
enough accumulated evidence that the KLN theorem works also at finite temperature 



H, |4|,|7|. 



Hence, the resummation of Braaten and Pisarski seems really to work well and 
certainly oversteps its original goals. Even in the case when the weak coupling limit 
is not justified, this separation between soft and hard scales can be used as a math- 
ematical trick, which is certainly more elegant than any other method. I only want 
for proof the numerous computations using these techniques [32-36,42,43,58,65,69]. 

5 Conclusion 

The list of references shows by its length that the TFT is now a mature field, with 
a growing activity over the last few years. A number of important questions have 
already been answered, such as for instance the comparisons between the different 
formalisms. The main difficulties of the TFTs are of course infrared singularities. A 
lot of effort has been done and the situation can be grossly summarized as follows: 
one can make reliable perturbative calculations when probing the thermal system at 
a scale T. By using the resummation of Braaten and Pisarski one can go down to the 
scale gT. The next step is to probe the magnetic mass scale, g'^T, onto which very 
little is known. 

There are many issues which are still unresolved. For instance, nobody knows 
exactly what the behavior of the running coupling constant is at finite temperature, 
and even if the question is relevant or not [ff^ . 

Finally, I could not end this paper without mentioning the interesting connection 
between thermal and topological field theories as, for instance, the fact that the 
Minkowski vacuum agrees with a thermal state for an accelerated observer ||75|. 
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Figure Captions 



Fig. 


1 


The real-time contour in the complex-time plane. 


Fig. 


2 


The deformed contour for calculating (3.2) in the ITF. 


Fig. 


3 


The two topologies for the self-energy at 2-loop order. 


Fig. 


4 


A "daisy" diagram with an arbitrary number of loops attached to the tadpole. 


Fig. 


5 


The effective propagators and vertices. 


Fig. 


6 


The dispersion relations for the bosonic modes. 


Fig. 


7 


The dispersion relations for the fermionic modes. 


Fig. 


8 


The resummed diagrams for the soft gluon damping rates. 
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